Non-universal ordering of spin and charge in stripe phases 
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We study the interplay of topological excitations in stripe phases: charge dislocations, charge 
loops, and spin vortices. In two dimensions these defects interact logarithmically on large distances. 
Using a renormalization-group analysis in the Coulomb gas representation of these defects, we calcu- 
late the phase diagram and the critical properties of the transitions. Depending on the interaction 
parameters, spin and charge order can disappear at a single transition or in a sequence of two 
transitions (spin-charge separation). These transitions are non-universal with continuously varying 
critical exponents. We also determine the nature of the points where three phases coexist. 
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High-T c compounds fascinate not only because of su- 
perconductivity but also because of a variety of concur- 
ring orders. In particular, theoretical Q and experimen- 
tal Q evidence has been found for stripes. Holes which 
are induced by doping condense into arrays of parallel 
rivers in the Cu02 layers. Within a layer, each river acts 
as a boundary between antiferromagnetic domains with 
opposite sublattice magnetization. Thus, stripes are a 
combined charge- and spin-density wave. 

Based on charge density and magnetization as order 
parameters, it is instructive to analyze the interplay of or- 
ders in the framework of a Landau theory H . However, in 
low dimensional structures fluctuations can be crucial for 
the nature of phases and of phase transitions. In particu- 
lar, fluctuations play a central role (i) for spin-charge sep- 
aration, i.e., the phenomenon that charge order emerges 
at higher temperatures than spin order (as observed in 
cuprates Q] as well as in nickelates d)) an d (ii) for the 
anomalous properties of the cuprates near optimum dop- 
ing H . To account for collective low-energy excitations of 
the electronic system, continuous deformations of perfect 
stripe order (spin waves or smooth stripe displacements) 
as well as topological defects (such as dislocations, vor- 
tices, or skyrmions) must be considered. The latter were 
found to induce transitions between various liquid-crystal 
like electronic phases |7). Besides transitions which are 
related to a degradation of the charge and spin struc- 
ture factors, Zaanen et al. |8| have suggested a further 
transition characterized by a less accessible, intrinsically 
topological order. 

In this Letter we present a paradigmatic model which 
is amenable to a largely analytical analysis of the in- 
terplay between charge and spin orders. Motivated by 
the weakly coupled layered structure of the materials, we 
restrict our analysis to two dimensions. Since the rel- 
evant materials typically have a planar spin anisotropy, 
the out-of-plane component of the spins is neglected. As- 
suming quantum fluctuations to be weak in comparison 
to thermal fluctuations we treat all degrees of freedom as 
classical. In the framework of a renormalization-group 
approach we establish the phase diagram and the nature 



of the phase transitions. The latter are driven by three 
classes of topological defects (cf. Fig. [j]): mixed defects 
combining a dislocation in the charge-density wave with 
a half-vortex, charge loops (with a Burger's vector of two 
stripe spacings) and entire vortices 




FIG. 1: Different types of topological defects in the charge 
density (lines) and sublattice magnetization (bold arrows), 
(a) Combination of a stripe dislocation (dashed arrows cor- 
respond with the Burger's vectors) and a half-vortex, (b) a 
charge loop, and (c) a vortex. 

We identify four different phases (cf. Fig. ||), depend- 
ing on which types of topological defects proliferate. In 
phase I there are no free defects. In phase II only vor- 
tices, in phase III only charge loops and in phase IV 
all types of defects are present. We characterize these 
phases by the range of charge order (CO), spin order (SO) 
and collinear order (LO). The phases are separated by 
transition lines along which the correlation functions de- 
cay with non-universal, continuously varying exponents. 
Slightly above the temperature T\ of the transition out of 
phase I, the short-ranged orders have a correlation length 
£ ~ exp{c(T — T\)~ v } with v = i except for the triple 
points Pi j2 (cf. Fig. |]) with v = ~ . 

To be specific, we assume that the stripes are paral- 
lel to the y direction with a spacing a. In the ground 
state, the charge density is modulated with a wave 
vector q = (27r/a,0). In the absence of dislocations, 
the stripe conformations can be described by a single- 
valued displacement field u = (u, 0) and the domains be- 
tween the stripes can be labeled by the function 6(r) = 
q • (r - u) + Y,m^Q i^I e4mq ' (r ~ u) which increases by 27r 
across a stripe. The charge density is located at the do- 
main boundary and can be expressed as p(r) = j^d x 6{r) 
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FIG. 2: Phase diagram. Thin dashed lines show where de- 
fects would become relevant for vanishing fugacities. For finite 
bare fugacities M the transitions are located at the bold lines. 
Phase I has no free defects and therefore CO, SO, and LO ex- 
ist. In phase II, vortices proliferate and destroy SO. In phase 
III, free loops are present, destroying CO and SO. Eventually, 
all orders are absent in phase IV. 



where A denotes the charge per unit length per stripe. 
For XY spins cr(r) = {cos $(r), sin $(r)}, the magnetiza- 
tion (including the modulation by the antifcrromagnetic 
order within the domains and by the anti-phase bound- 
ary condition on the stripes) is captured by the angle 
$(r) = Q-r + #(r)/2 + 0(r). Here Q is the antiferromag- 
netic wave vector and <p describes the smoothly varying 
spin Goldstone modes. 

In the absence of topological defects low-energy exci- 
tations are wave-like and governed by the Hamiltonian 



J d 2 r[ J s (V0) 2 + J c (2ir/a) 2 (Vu) 2 } (1) 



with spin and charge stiffness constants J a (with a — s, 
c, respectively). The structure factors S p (r — r') = 
(5p(r)Sp(r')) and S a (r — r') = (<r(r) • cr(r')) decay al- 
gebraically, S p (r) ~ cos(q • r)r~ T,p with r\ p = 1/(2ttK c ) 
and 5 CT (r) - cos[(Q + f ) • r]^" with r\ a = 1/(2ttK s ) + 
1/(8ttK c ). We define the reduced stiffness constants 
K a = JJT. 

In analogy to the Kosterlitz-Thouless (KT) transi- 
tion |ic], [ll[], the presence of topological defects effects 
a screening of the stiffness constants. If neutral defect 
clusters unbind, these constants are renormalized to zero. 
For suitable values of the stiffness constants charge-loop 
pairs unbind as only type of defects (phase III). Then 
both S p and S a decay exponentially. To distinguish this 
phase from phase IV where all defects are free, an addi- 
tional correlation function is necessary. For this purpose, 
Zaanen et al. |8|, [l^] have suggested a highly nonlocal cor- 
relation function involving spin and charge. As a simpler 



and probably more physical alternative, we propose 

C\ ] {v-T , ) = 2{{a T -a T ,f)-l (2) 

which measures the spin collinearity and has the advan- 
tage of being local and being defined entirely in terms of 
spin variables. Collinearity is insensitive to continuous 
stripe displacements and to loops since the spin ground 
state remains collinear in the presence of such excitations. 
In the absence of topological defects in the spin sector (in 
phase I and III), spin waves lead to an algebraic decay 
of collinear order (LO), Cm (r) ~ cos(2Q • r)r~ ?, ii with 
77,1 = 2/(nK s ). 

We now calculate the mutual screening of the defects 
which, in two dimensions, interact logarithmically on 
large distances. Therefore, their interactions can be cap- 
tured by a Coulomb-gas formulation 



H top = ^2 ^2 ^Ja"fn ia m ia In — 

i^j a— s,c 



(3) 



for topological vector charges (mi S ,mi c ) at positions rj. 
Additional core energies of such vector charges give rise 
to bare fugacities yt = exp{—TrC(K s m 2 s + K c m 2 c )}. C 
is a constant of order unity ||. Vortices have a vec- 
tor charge (m B ,m c ) — (±1,0) and a fugacity y v = 
exp(— TtCK s ), and loops are characterized by (m s ,m c ) — 
(0, ±2) and y\ = exp(—4TtCK c ). Dislocations have vec- 
tor charges |m c | = 1 and \m s \ = ^ and fugacities 
y t \ = cxp{— nC(K c + K s /4)}. Defects of higher charges 
are negligible for the critical aspects of the phase dia- 
gram. 

We follow the renormalization approach developed by 
Kosterlitz |ll]] for the usual scalar Coulomb gas and gen- 
eralized to vector gases in different contexts H& The 
charges are assigned a hard-core cutoff a which is in- 
creased to a = ae dl under an infinitesimal coarse grain- 
ing. Thereby, pairs of vector charges at a distance a 
annihilate each other if they have opposite charges, oth- 
erwise they recombine to a single non-vanishing vector 
charge. This coarse graining procedure leads to scale de- 
pendent stiffness constants and fugacities, for which we 
obtain the flow equations 



dK: 



dl 

dKr 1 



= 2ir 3 (2y 2 v + y 2 ), 
= 87r 3 (2y 2 + y 2 ), 



dl 

^ = (2-7^)^ + 27^ 



d! 



^ = (2 - A7rK c ) yi + 2iry 2 



d) 



(4a) 
(4b) 
(4c) 
(4d) 



dl 



(2 - ^(K s + 4kS) Vd + 27r(y v + Stt)w.(4e) 



They are invariant under the mapping 
(AK c ,K s ,y h y v ) -> (K B , 4K C , y v , y{). A qualitative 
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understanding of the phase diagram can be obtained 
from the rescaling contributions linear in the fugacities to 
the flow equations. Vortices proliferate for irK s < 2 and 
loops for AttK c < 2. Furthermore, for j(K s + AK C ) < 2 
dislocations become relevant. The corresponding borders 
are shown in figure Fig. || as dashed lines. 

In order to include the terms quadratic in the fugac- 
ities, the flow equations have to be integrated numeri- 
cally. In phase I all fugacities tend to zero as I — > oo and 
K s and K c are renormalized to finite values K E (oo) and 
K c (oo) smaller than the initial ones. The defects can be 
considered as being bound in local clusters of vanishing 
total charge. Their fluctuations enhance the wave exci- 
tations such that CO, SO, and LO are quasi- long ranged: 
S p , S fj , and C\\ decay algebraically like in the absence of 
defects. The exponents rj are given by the expressions 
given after Eq. ([!]) if the bare stiffness constants are re- 
placed by the renormalized ones. Since the boundary of 
phase I flows to a line of fixed points (dashed lines in 
Fig. H for zero fugacities) with continuously varying val- 
ues of the renormalized stiffness constants, the transition 
in non-universal with exponents given in Tab. S. 
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TABLE I: Values of the exponents 77, v at the phase transi- 
tions. Orders that become short-ranged at a transition line 
(77 jumps to infinity) are marked by an asterisk. 

Outside phase I, at least one fugacity increases with the 
scale. Since our flow equations are valid only for small fu- 
gacities, they can be evaluated outside the ordered phase 
only up to a finite scale I* where some fugacity becomes 
of the order of unity. The divergence of a fugacity drives 
one or both stiffness constants to zero. A divergence of 
the renormalized exponent 77 then signals short-ranged 
order with a finite correlation length, which scales like 
£ ~ ae l sufficiently close to phase I. Although the diver- 
gence of fugacities hints at the nature of the phases II, III, 
and IV, the precise shape of these phases is determined 
by strong coupling (large fugacity) regimes beyond the 
validity of Eqs. (f|). 

In phase II, vortices proliferate and the spin stiffness 
K s is renormalized to zero. This leads to a destruction of 
the spin and the collinear order, since the exponents r] a 
and 7711 are infinite. The corresponding correlation func- 
tions decay exponentially, S a ~ C\\ ~ cxp(— r/£). Since 
K s is renormalized to zero, the interaction in the spin 
sector effectively breaks down. Nevertheless, dislocations 



and loops are still coupled in the charge sector and can 
remain bound for sufficiently large K c Then also 

charge order remains quasi-long ranged with finite r\ p . 

In phase III loop pairs unbind and the charge stiffness 
K c is renormalized to zero, leading to r\ a = r\ p = 00 and 
therefore to a short ranged spin and charge order. The 
spin order is destroyed because of the arbitrarily large 
fluctuations of the domain walls. Since K s is screened 
only by bound vortices and dislocations, ryy remains finite 
and quasi-long ranged collinear order is preserved. 

In phase IV the proliferating dislocations obviously 
render all correlations short ranged, S a ~ S p ~ Cn ~ 
exp(-r/£). 

As pointed out above the precise location of the tran- 
sitions between phases II and IV and phases III and IV 
cannot be obtained by the weak coupling (low fugacity) 
flow equations. Nevertheless, we now present qualitative 
arguments for their location, focusing on the example of 
the transition II/IV. In the limit K s (0) = the inter- 
action in the charge sector is simply switched off. Nev- 
ertheless, loops and dislocations interact in the charge 
sector; the former stronger than the latter. Therefore, 
the transition (point Q\ with K c (po) = 2/tt) is driven by 
the unbinding of the dislocations. At this point loops are 
irrelevant and the transition belongs to the KT univer- 
sality class. When the bare K s is increased, it is renor- 
malized to zero by free vortices and the transition occurs 
practically at the same value of K c as long as K s < 2/ir. 
When K s increases to the value where vortices start to 
bind, the transition line II/IV reaches the point P\ at 
a smaller value of K c since dislocations start to be sta- 
bilized by the additional interaction in the spin sector. 
The situation for the transition line between the phases 
III and IV is similar, in this case the interaction of the 
charge components of the mixed defects is screened by 
the proliferating charge loops. 

At the triple points P\,i the critical properties are gov- 
erned by the balanced competition between two types of 
defects which can recombine into each other. Therefore, 
these points do not belong to the universality class of the 
KT transition. Due to the symmetry of the flow equa- 
tions we concentrate on point P2 . There one can neglect 
y v to a good approximation and the conditions y\ = y& 
and K s — 12K C (thin dotted line in Fig. || top) are con- 
served under the flow. In this case the flow equations 
for Y ;= 2-7r\/6yd and the small variable x := —k, 2 

read = Y 2 and ^ = xY + -^Y 2 . Following Refs. 
| , these differential equations can be solved analyti- 
cally and result in a critical exponent v — |. 

Thus, the exponent v is discontinuous along the border 
of phase I. While this discontinuity describes the diver- 
gence of £ in the thermodynamic limit, it is unlikely to 
be seen in samples of finite size. For large but finite sys- 
tems sizes, £ appears to diverge with an effective expo- 
nent i/ e fi which appears to vary continuously in the range 
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0.4 < i/ e g < 0.5. We have analyzed the divergence of 
£ from a numerical integration of the flow equations up 
to Z max = f 3 which corresponds to astronomically large 
scales. The resulting v c g (cf. Fig. ||) shows a notch of 
a relatively small width which can become substantially 
larger for smaller values of Z ma x- 



turn critical points at zero temperature and the influence 



of disorder (see also 16 ) on the nature of ordering. 
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FIG. 3: Top: magnification of the phase diagram near point 
P-2- In this region, phase III is extremely narrow and masked 
by the transition line to the right of P2 . Bottom: Numerically 
calculated effective exponent v e s for i max = 10 3 along the 
boundary of phase I. 

From the phase diagram of our model several scenarios 
of spin-charge separation are possible. At low tempera- 
tures the stripe phase I is realized, (i) For J c > |j s , 
there occur two distinct transitions with increasing tem- 
perature: first into phase II (loss of SO and LO), then 
into phase IV (loss of CO). This scenario is observed in 
experiments ||, ||. (ii) For -^J s < J c < |J S there is a 
single transition from phase I to phase IV where all orders 
disappear simultaneously, (iii) For J c < yjJs, CO and 
SO disappear at the same temperature of the phase tran- 
sition I /III while LO disappears only at an even higher 
temperature at the transition III/IV. The relation be- 
tween J s and J c - which determines the scenario - can 
be tuned by doping. With increasing doping, the shrink- 
ing stripe distance should lead to a significant increase of 
J c , whereas J s should change only weakly since it is es- 
sentially determined by the antiferromagnetic exchange 
coupling. Therefore, in principle, the last scenario could 
be realized and the transition III/IV should be detectable 
for example by polarized neutron scattering. 

In conclusion, we have examined a model for coupled 
spin and charge order in stripe phases. We have shown 
that the transition between the various phases are non- 
universal. We have identified collinearity as a physical 
quantity that allows to discriminate between phases III 
and IV. Several interesting questions invite to future in- 
vestigations, in particular the properties of related quan- 
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